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Abstract 

Quantum mechanical tunneling across smooth double barrier potentials modeled using Gaussian 
functions, is analyzed numerically and by using the WKB approximation. The transmission prob- 
ability, resonances as a function of incident particle energy, and their dependence on the barrier 
parameters are obtained for various cases. We also discuss the tunneling time, for which we obtain 
generalizations of the known results for rectangular barriers. 
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I. INTRODUCTION 



The usual example of quantum mechanical tunneling is the rectangular barrier in one 
dimension.- - - Curious students might wonder what happens if we consider smooth barriers. 
Do the tunneling results remain the same? Are there quantitative or qualitative differences? 
Students might also wonder what happens if there is more than one barrier. To answer 
these questions (particularly for smooth barriers), numerical methods and approximations 
are essential because little can be done analytically. In this paper we discuss tunneling across 
smooth double barriers of various types. 

Double barrier potentials arise in many diverse areas. We give a few examples in the 
following. 

Quantum heterostructures^ Semiconductor heterostructures are layered, thin (about 100 
nanometers or less) sandwich structures made with different semiconductor materials (for 
example, GaAs between two AlAs layers). The existence of junctions of different materials 
is the reason for the occurrence of sequences of wells and barriers in such structures. Semi- 
conductor heterostructures have made it possible to control the parameters inside crystals 
and devices such as band gaps and the effective masses of charge carriers and their mobili- 
ties, refractive indices, and electron energy spectrum. Heterostructure electronics are widely 
used in many areas, including laser-based telecommunication systems, light-emitting diodes, 
bipolar transistors, and low-noise high-electron-mobility transistors for high- frequency ap- 
plications including satellite television.- 

High energy physics. A barrier penetration model has been developed for heavy ion 
fusion, 6 taking into account the realistic features of the Coulomb potential for the case of 
nuclei. The typical barrier shapes encountered for nuclei conform to double barriers. 

Extra dimensions. Another area where such potentials arise is in the physics of extra 
dimensions. In particular, they appear in the context of localization of particles and fields 
on a four- dimensional hypersurface, known as the brane. In these brane-world models 7 
localized particles in four dimensions are viewed as bound or quasi-bound (resonant) states 
in a double barrier effective potential spread across the extra dimensional coordinate. 

Nonlinear Schrddinger equation. The nonlinear Schrodinger equation has been used in 
investigations of Bose-Einstein condensates to probe macroscopic quantum tunneling^ - — 
and gravity surface waves in fluids, among many other applications. The methods we will 
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discuss here may be extended, with minor modifications, to the Gross-Pitaevskii equation, 
a nonlinear Schrodinger equation arising in the study of many body systems. Above barrier 
reflection and tunneling in the nonlinear Schrodinger equation is discussed in Ref. si The 
double Gaussian barrier, as well as the rectangular barrier, has been used as model potentials 
in studies^ 1 ^ of the nonlinear Schrodinger equation. However, the definition of the tunneling 
coefficient is ambiguous for weak nonlinearity, because the principle of superposition does 
not strictly hold. 

Apart from the applications we have discussed, we mention some recent pedagogical 
articles on double barrier tunneling. Numerical solutions for quantum tunneling through 
rectangular double barrier heterostructures have been explored.™ Further, a detailed analysis 
of the propagation of wave packets across double barriers has been performed.— Three- 
dimensional s-wave tunneling for double barrier potentials has also been investigated.— 

The paper is organized as follows. In Sec. II we introduce the various types of double 
barriers that will be discussed. Section III discusses the theory of double barrier tunneling, 
and the various methods we will use. Numerical and semiclassical analyses of the smooth 
barriers mentioned in Sec. II are carried out in Sec. IV. In Sec. V we determine the tunneling 
time across smooth double barriers. Section VI summarizes our results and mentions some 
possibilities for future studies. 



II. TUNNELING ACROSS DOUBLE BARRIERS 



A simple way to model double barrier potentials is to add a term shifted in position to 
the single barrier potential function, that is, 

V{X) = ^singleO, Vi, Wi) + V sing i e (x - a, V 2 , W 2 ), (1) 

where Vi j2 and u>i j2 are the height and width of the two barriers. The separation between 
the barrier heights is a. 

The simplest and the easiest barrier to analyze is the rectangular double potential barrier, 
whose parameters are defined in Fig. [TJ We have 



V sinK i e (x,V 1 ,w 1 ) 



x < 

Vi 0<X<W! 

x > W\ 



(2) 



FIG. 1: Rectangular double barrier (solid line) of widths w\ and w 2 separated by the distance a. 
The heights of the two single barrier potentials are V± and V 2 . Also shown is a double Gaussian 
barrier (dashed line). 
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Similarly, the Gaussian double barrier (see Fig. [T]) is written using 



V single (x } V u a l ) = V 1 e- x2 / 2 ^ 



(3) 



and hence V(x) = V singlc (x, V 1 , at) + V sing y e (x - a, V 2 , a 2 ). 

The time- independent Schrodinger equation for a one- dimensional potential V(x) is 

^ + f^-V^ = 0, (4) 

where m is the mass of the particle, E is its energy, and if) is the wavefunction. The 
Schrodinger equation for a rectangular double barrier can be readily solved for if), giving 
five solutions in the five regions. The solutions and their derivatives can be matched at the 
boundaries between two regions to obtain a relation between the incoming and outgoing 
waves. The transmission amplitude isM 

T = 4{e ikl( - Wl+W2) (2cosh(k 3 w 2 ) - ifci_3sinh(fc 3 W2))(2cosh(fc 2 wi)) - iJfei_ 2 sinh(Jfe 2 itfi) 

+ e ikl{9+b - Wl) k 1+3 k 1+2 smh(k 2 w l ) sinh(A;3w; 2 )]- 1 , (5) 
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where W\ and w 2 are given in Eq. ([I]),©. The other parameters are defined as 



g = wi+w 2 , b = wi + a + w 2 (6a) 

k x = V2mE/h, k 2 = V 2m (^i - E)/h (6b) 

k 3 = y/2m(V 2 -E)/h), k i±j =-^±-t : (i,j = 1,2,3). (6c) 

rCj r\i<i 

The transmission coefficient (probability) is T = T*T. 



For the rectangular double barrier, for which we have analytical solutions, we match ijj and 
dijj/dx at the discontinuities to obtain the relations between the wavefunction amplitudes, 
that is, A and C. For a smooth barrier, the ratio of A and C can be obtained by solving 
Eq. (|4j) numerically. The ratio of the maxima on the right-hand side of the barrier to the 
incident amplitude on the left-hand side of the barrier can be found from the plot of ip versus 
x. 

The analysis can be extended to the arbitrary potential functions.- 1 ^ However, exact 
analytical solutions of Eq. (JJ]) are usually difficult or impossible to find. In particular, it is 
impossible to exactly solve the Schrodinger equation for the Gaussian potential. To obtain 
useful results, we either find numerical solutions or use the WKB approximation. 

To solve the Schrodinger equation numerically, we write Eq. (jl]) as a system of two first- 
order linear ordinary differential equations with <f>i — ijj and <ft 2 = dip/dx. With these 
substitutions, Eq. fl4| becomes 

0i = 02, (7a) 
2tt? 

<t> 2 = w (V(x)-E)<t> 1 . (7b) 

Equation (JTj) can be solved for (f>i to obtain ip(x). We used the fourth order Runge-Kutta 
algorithm. 

We compute the wavefunction on both sides of the barrier for different energies in the 
required range (up to twice the maximum potential height of the barriers). Because the 
potentials are rapidly decaying and we are mainly interested in the barrier penetration 
properties, the use of asymptotic solutions for ip to calculate T is justified. To the far right 
and far left of the well, the forms of ip are, 

bright -> Ce lkx , Vieft Ae lkx + Be~ lkx . (8) 
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The transmission coefficient is calculated from 



r 



\A\< 



(9) 



The Wentzel-Kramers-Brillouin (WKB) approximation^^ can be used to obtain partially 
analytic results for an arbitrary potential well. In this approximation the wavefunction if) is 
expressed as a power series in h. We write the wavefunction in the form 



if) = A(x)e 



(10) 



If we assume the potential is slowly varying and neglect the second derivatives of and A, 
the solutions of the Schrodinger equation can be expressed as 



iP = A\p\ 



-l/2 e ±if\p\dx/h 



V 



yj2m\E -V{x)\. 



(11) 



The limits of integration in Eq. ( TTTj) are determined by the classical turning points at which 
E = V(x). At these points the semiclassical approximation fails because the wavefunction 
becomes infinite. To the far left and to the far right of the barriers, the WKB method gives 
reasonably good approximations for if)(x). The matching of if) near the classical turning 
points is done by the use of special functions. The asymptotic forms of these functions result 
in a net phase shift of 7r/4 when going from a classically forbidden region to a classically 
allowed region or vice-versa.- 1 ^ 



E / \ 




J \ 





FIG. 2: An arbitrary double barrier potential V(x) with classical turning points at x\, X2, x^, and 
X4. The regions x\ < x < X2 and X3 < x < X4 are classically forbidden, but the particle can tunnel 
through these regions quantum mechanically. 
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For a general single barrier potential, the transmission probability is 

-2 



16 



r 



1 T 

T + I 



T 2 /(l + T 2 /4) 2 



where 



T = exp 



X2 

y/2m(V(x) - E)dx/h 



■t'i 



(12) 



(13) 



and X\ and £2 are the classical turning points. The WKB approximation introduces errors 
which are consistent with T being small for the single barrier case. Thus T — T 2 , or 



T — exp 



X 2 

y/2m(V(x) - E) dx 



(14) 



For a double barrier we obtain four classical turning points instead of two as evident from 
Fig. [2J Starting with the wave function on the far right of the barrier (region V), we have 



ipv = Ap 1 exp ( i J p/Tidx + i— 



7T 



.X'4 



(15) 



Let T\, T 2 , and T 3 be defined as 



T\ = exp 
T 3 = exp(- 



X4 



\p\/hdx 



p/hdx 



X2 



\p\/hdx) 



(16a) 
(16b) 
(16c) 



We expand the complex exponentials using Euler's identity, and then match the real and 
imaginary parts of if) near x < £4 and x > £4, and then near x < £3 and x > £3 (by applying 
the connecting formulae of the WKB method) and obtain 



in 



iAp 



-1/2 



Ti_l_ 
4 T x 



exp ( i / p/hdx + 



Similarly, we find 



A 



+ 



'Ti 1 
+ ( T + -]exp 



C3 + C4 iT 3 C±-iT 3 C 3 
*T 3 4 

C 3 + C 4 iT 3 C 4 - iT 3 C 3 



p/hdx — 



(17) 



in 



4 



exp 1 2 



exp 



j-i 



p/hdx + 
p/hdx 



(18) 
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The second term can be identified with the wave incident from the left, and yields the 
transmission coefficient 



T 



Vv 


2 


T^{C^ — C3) 


V'incident 




4 



+ 



where 



C 3 = (l/T 1 -T 1 /4)exp(iT a ) 
C 4 = (Tx/4 + 1/TO exp(-iT 2 



(19) 



(20a) 
(20b) 



and = y/2m[E — V"(x)]. 

Equation (Tl9|) can be applied to a Gaussian double barrier with the assumptions that 



E < Vx, V 2 , a ^> <7i, cr 2 , and a — 3ai — 3a 2 > 0. Hence, £i )2 = ±ct"i a/2 (log Vi — log-E 1 ), and 
X3 t 4 = a ± o" 2 A/2(log V 2 — log-E). The integrations in the WKB method were performed 
using Simpson's rule. 



III. NUMERICAL AND SEMICLASSICAL ANALYSIS FOR MODEL SMOOTH 
BARRIERS 



A. Numerical Analysis 

The T versus E plots in Fig. [3] show the variation of the transmission probability with 
the energy of an electron incident on a one-dimensional Gaussian double barrier. Figure [3] 
illustrates how the transmission coefficient changes due to the presence of two barriers, a 
change in the barrier separation, and a change in the height of one barrier. The energy of the 
electron is assumed to be of the order of a few electron volts. In Fig. [3]^a) the two barriers are 
very close to each other and the number of resonances is less (Fig. [3^b)) than obtained for 
larger a. In Fig. [3](c) the separation between the barriers is increased,the well becomes wider 
than the well in Fig. [3](a), and consequently, the number of resonances increases (Fig. [3](d)). 
Figure [3]^e) shows two barriers that are so close that they appear to be a single continuous 
barrier, and hence no resonances (Fig. [3](f)) are seen. In Fig. [3](g) we consider asymmetrical 
barriers for which V\ 7^ V 2 and/or o\ ^ a 2 with a\V\ = o^i- We note that the effect is 
predominantly that of a single barrier though resonances do appear (Fig. [3](h)), signaling 
the presence of a second barrier of much smaller height. 



S 







Energy in eV 



FIG. 3: Plot of the potential (shaded) and wave function. The plots of the wave function are for an 
incident energy of 4 eV. (a) V\ = V2 = 4cV, g\ = <r 2 = 0.2 nm, a = lnm. (b) Transmission probability 
for the barrier in (a), (c) V\ = V2 = 4eV, u\ = <r 2 = 0.2 nm, a increased to 3nm. (d) Transmission 
probability for the barrier in (c) - the number of resonances is greater due to the larger separation-to-width 
ratio, (e) V\ = V2 = 4eV, u\ = cr 2 = 0.2 nm, a reduced to 0.4 nm so that the two Gaussian potentials 
almost merge into a single barrier, (f) Transmission probability for the barrier in (e)- the plot resembles 
the transmission probability of a single Gaussian barrier, (g) V\ = 1, V2 = 4cV, <j\ = 0.4, a 2 = 0.1 nm, 
a = 3 nm. (h) Transmission probability for the barrier in (g)~ asymmetrical barriers have a marked effect 
on the transmission probability. 




Energy) in eV ) Energy (in eV) 



FIG. 4: Comparison of transmission probabilities for the rectangular and the Gaussian double barrier, (a) 
Square : V\ = 4, V 2 = 4 eV, wi = 0.6, w 2 = 0.6, a = 0.75 nm (b) Gaussian : Vi = 4, V 2 = 4eV,ai = 0.6, cr 2 = 
0.6, a = 4.35 nm (c) Square : Vi=A,V 2 = A eV, w 1 = 0.2, w 2 = 0.8, a = 1 nm (d) Gaussian : Vi = 4, V 2 = 4 
eV, tTi = 0.2, U2 = 0.8, a — 4 nm (c) Square : Vi = 4, V 2 = 8 eV, wi = 0.2, w 2 — 0.1, a — 4 nm (f) Gaussian 
: Vi = 4, V 2 = 8 eV, (Ti = 0.2, o 2 = 0.1, a = 4.9 nm The transmission probability asymptotically approaches 
unity for the Gaussian barrier due to the smoothness of the Gaussian barrier potential. 

Next, we compare the transmission coefficients for the Gaussian and the rectangular 
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double barrier (see Fig.|4j. We choose a(Gaussian) = 3wi + 3w2 + a(rect), o\ = w±, and o<i = 
u>2- We also vary the height of the barriers keeping the (width x height) for both the barriers 
the same. Figure H^a),(c),(e) shows the transmission coefficient for a double rectangular 
barrier, and Fig. H]^b),(d),(f) shows the same for a double Gaussian. The qualitative nature 
of the plots for the rectangular and the Gaussian barriers is similar. Tunneling occurs for 
both barriers, though there are differences that may be used to determine the barrier shape 
from plots of T versus E. 

The number of peaks in T is large for both cases. However, the peaks and valleys are 
more prominent for the rectangular barrier. The transmission coefficient for the Gaussian 
does not reach unity (except asymptotically) and the minima are also at lower values. The 
cause of this quantitative difference is the smoothness of the Gaussian barrier. 

The following observations are valid if a/a is large (> 4). The T versus E plot (Fig. H]) 
depends strongly on the product of width and height for any barrier. If we keep the width 
x height constant, the graphs are qualitatively similar. 

The plots show a large number of valleys and peaks, corresponding to resonant states. 
Resonant states are a result of the destructive interference between waves reflected from the 
two barriers. The valleys in the T versus E plots correspond to quasi-bound states. As the 
separation-to-width ratio of the barriers is increased, the number of resonances increases, 
which means that the well can accommodate a larger number of quasi-bound states. If the 
well width is made indefinitely wide, the resonant states merge into a continuum. 

The peaks become sharper as the height of the barriers is increased; i.e. the full width at 
half maximum of the peaks decreases. If A is the full width at half maximum of the peak 
in T(E), then A decreases as Vq an d V\ increase. The peaks also become sharper if the 
widths of the barriers are increased, which makes it difficult to obtain numerical solutions. 
The lower the energy at which the peak occurs, the steeper is the peak. If E <C Vq, V±, A 
is very small and vice versa. 

B. WKB Analysis 

The WKB method provides a way of obtaining analytical expressions for the transmission 
coefficient. In Figure |5] we compare the results obtained for the transmission probability 
using the semiclassical analysis method and the numerical solutions Figures EJ^a) and (b) 
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differ because of the larger value of a in Fig. GO^b) (all other parameters are the same). In 
Fig. |5](c) the heights are halved, and the widths are doubled compared to Fig. [5(b). In 
Fig. 0(d) the Gaussian barrier is asymmetrical with cri/cr2 — 6. 

For £ C Vi the values of T are almost the same, or at least of the same order of 
magnitude using either the numerical method or the WKB approximation. The greatest 
deviations from the numerical solution are observed for energies near Vq. This deviation is 
expected because the WKB approximation is not valid near the classical turning points. The 
difference between the WKB method and the numerical solution of Schrodinger's equation 
also becomes prominent when the separation between the barriers becomes comparable to the 
width of the barriers. In the calculation of the turning points we neglected the contribution 
of the first barrier to the turning points of the second barrier and vice versa so that we 
could obtain analytic expressions for the classical turning points; otherwise, we would have 
to solve the transcendental equation V{x) = \V\ exp(— x 2 /2a\) + V 2 exp(— (x — a) 2 /2a 2)} = E 
for each value of E. Thus the difference between the semiclassical results and the numerical 
solutions for the transmission probability becomes prominent when the barriers are close to 
each other, but the assumption yields acceptable results in return for a shorter computation 
time. 

The WKB approximation becomes more accurate as the separation between the barriers 
increases. The numerical values of T(E) versus E in Fig. [5] do not match the WKB approx- 
imation results for energies E > min(Vi, V2). For asymmetrical barriers with non-identical 
potential heights, the WKB approximation breaks down for E > rnin(Vi, V2) because the 
formulation we have used involves four classical turning points. The potential function re- 
duces to a single barrier with two turning points for E > min(Vi, V2). However, the WKB 
results agree quite well for E < min(Vi, Vq). 

The advantages of the WKB method are that it requires about an order of magnitude less 
computation time than the direct numerical solution of Schrodinger's equation. The main 
limitation of the WKB method is that it relies heavily on the turning points, and hence is 
difficult to apply to energies greater than the potential height of any of the barriers. 
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0.5 1 1.5 2 2.5 3 1 2 3 4 



Energy in eV Energy in eV 

FIG. 5: Comparison of the transmission probabilities calculated by numerical solution of 
Schrodinger equation (gray) and the WKB approximation (black) for a Gaussian double bar- 
rier with parameters: (a) V\ = Vi = 4eV, o\ = 02 = 0.2 nm, a = 1 nm, (b) V\ = V2 = 4eV, 
u\ = (72 = 0.2 nm, a is increased to 4nm. (c) V\ = V2 = 2eV, o\ = 02 = 0.4 nm, a = 4nm. 
The potential barrier height is larger for this barrier, (d) V\ = V2 = 4eV, a\ = 0.6, 02 = 0.1 nm, 
a = 4nm. In this case we consider an asymmetrical barrier. The results obtained by the two 
methods are in good agreement. 
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IV. TUNNELING TIME 



With the advent of the fabrication of nanometer semi-conductor structures, the calcula- 
tion of tunneling time has acquired more importance.— For one-dimensional single barriers, 
the tunneling time becomes independent of the barrier width for opaque barriers. This 
"Hartman effect"— predicts superluminal and arbitrarily large group velocities inside suffi- 
ciently long barriers. There is some controversy concerning the definition of the tunneling 
time.— The expression for the tunneling time we have used is variously named the "group 
delay" or "phase time." 

If we assume a single rectangular barrier, with the potential 

{V Q (0 < x < L) 
(21) 
elsewhere, 

the wavefunction may be written as^ 

exp(ikx) + \/7Zexp(—ikx + ry) (x < 0) 

V>i(z) (0 < x < L) (22) 

V / exp(ikx + ia) (x > L). 



ip(x) 



where T and TZ are the transmission and reflection probabilities, k the wave vector, and a 
and 7 are phase constants dependent on k. If we apply the stationary phase approximation 
to the peak of the transmitted wave packet, d{axg(ip(x, t))}/dk = 0, we find 

da 1 dE „ . 

Tk ~ nTk* = °- (23) 

Thus, the temporal delay may be defined as 

From the transmission amplitude for a single barrier we find the phase shift to be a = 
— tan _1 ((/s; 2 — k 2 ) tanh(fcL)/2ft/c), where k = \/lvn{y§ — E)/% and k = \/2mE /Ti. By using 
Eq. (T51]) we have 

t » ^ & < 25 > 

where the subscript g denotes group delay. 
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For a rectangular double barrier the transmission amplitude is expressed by Eq. fl5]), and 
the phase time is^^i 

d 

r g = h—axg[Texp(ikib)]. (26) 

In the limit of k 2 Wi and /C3W2 — > 00, the phase time becomes independent of the barrier 
width and separation: 

T > = (27) 

except at resonance and anti-resonance. 

To illustrate what happens at resonance and anti-resonance,— we take E = Vi/2 and 
Vq — Vi — V 2 - When the phase shift 2k x a = 2mn {m = 1,2,3, . . .) (resonance), the group 
delay becomes 

res _ (1+jRo) « / 98 x 



where v = p/m = v2rnE /m, and Rq and T denote the reflection and transmission proba- 
bilities through a single barrier of width w±. At anti-resonance (2k\a = (2m + 1)tt), 

= (T^y^. (29) 

indicating that the phase time increases linearly with the inter-barrier separation, a, at both 
resonance and anti-resonance. 



A. Tunneling times for general double barrier potentials 

The wave function for x > 24 (Fig. [2]) is given in the WKB approximation by 

^ = A\p\-^ exp jf |p| rfx + i^j , (30) 
and the incident wavefunction is 
incident = A\p\~ 1 ' 2 {i(C 3 + C 4 )/T 3 + zT 3 (C 4 - C 3 )/4} exp 1 p/ft dx - m/^j . (31) 



In analogy with Eq. f )26l) . the group delay is found to be 



d 



X4 



(32) 



(V'V'Mncident) exp |p|/ftdz 

If we substitute and ^incident, into Eq. [32] we obtain 

7, -+- <7,1 _1 

(33) 



^3(^4 — C3) C4 + C3 
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and after some algebra 



d , 2tanT 2 

Tg= 8E tan T 2 ^ 



2 ^ 2rr «9T 2 tanT 2 9TA 
= TfT4ta^ Tl i SeC Ta flE " 2 ^T9E J • (34b) 

Equations [27H291 which are applicable for a rectangular barrier, can be recovered 
from Eq. (I34"l) as follows. We have T\ = exp(— /c 2 u>i) and T 2 = k\a. We let = 
tan~ 1 (2tan(T 2 )/T 1 2 ) and write r g = hd<p/dE. We find that the phase tunneling time through 
a rectangular double barrier of width w± and separation a is 

T g = 2 cos 2 0— exp(2fc 2 u>i)(a sec 2 (fcia)/A; 1 — 2w\ tan(A; 1 a)/A; 2 ). (35) 

To illustrate the behavior of r g near a resonance, we assume i? = V\/2 so that &i = A; 2 . At 
resonance, 2/^a = 2m7r, where m is an integer. The group delay at resonance assumes the 
form 

r res = 2e 2k2Wl a/v w (1 + R )a/T v, (36) 

because cos0 = 1 and tanfcia = 0. Here, i?o = (V^i — ^i/4) 2 /(l/Ti + Ti/4) 2 and T = 
(1/Ti + Ti/4)~ 2 are the reflection and transmission probabilities through a single barrier 
((1 + R )/T = 2{T{ 2 + T 2 /16) w 2e fc2W1 )- 

Similarly at anti-resonance (2/cia = (2m + l)vr), m an integer, we have 

T ^ = To a/v(l + R ). (37) 

We next obtain r g for a Gaussian double barrier. The relevant integrations for obtain- 
ing T 3 are performed numerically. Figure [6] shows the variation of r g with a for fixed a. 
Figures [7(a) and (b) show how r g varies with a for fixed a. 

We see that the group delay increases with an increase in the barrier width up to a certain 
value and then falls off rapidly on widening the barrier further. The peak in the group delay 
versus barrier width shifts toward the greater widths at higher energies. The maxima in r g 
shift to greater values of a at higher energies. The group delay decreases with increasing 
Vo, and for fixed Vq r g increases with E if other parameters remain constant. The tunneling 
time is more or less independent of the interbarrier distance as is evident from Fig. [7J which 
is consistent with the Hartman effect. 
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FIG. 6: Variation of tunneling time with barrier width a for a double Gaussian barrier, with (a) 
Vq = 4eV; a = 8nm (b) Vq = 12 eV; a = 8nm, for various energies. The tunneling time increases 
up to a certain value of a and then decreases rapidly. The maximum in the curve shifts to larger 
values of the barrier width for higher energies. 



(b) 



o 

a 

D 

o 




FIG. 7: Variation of the tunneling time through a double Gaussian barrier with inter-barrier 
distance a for (a) Vq = 6eV, a\ = 0.2 nm, and (b) Vq = 1.50 eV, o\ = 0.8 nm. Note that the 
tunneling time varies over several orders of magnitude for various energies of the incident particle. 
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V. DOUBLE BARRIERS IN HETEROSTRUCTURES 



Double barriers of either type (smooth or square) appear in various physical situations 
(as outlined in Sec. I). We now discuss in some detail, how they may arise in the context of 
one such scenario-semiconductor heterostructures. 

FIG. 8: Semiconductor heterostructures and superlattices: (i) Energy band diagram of an n- 
N semiconductor heterojunction; (ii) Hetero junctions can be of three types classified accord- 
ing to the lineup of the conduction and valence bands: straddled(£'c'2 > Eci,Ev2 < Ey x ), 
staggered(-E(72 > Ecx,Eci > Ey2 > Eyi), and broken(£'c2, Ey2 > Eci)) (in) An unbiased 
GaAs/AlGaAs superlattice and the corresponding potential energy diagram; (iv) A biased super- 
lattice created by the application of an electric field, which creates a gradient in the potential 
function. 
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(iii) 

A heterojunction (Fig. |8]d),(ii)) is the interface between two dissimilar solid state ma- 
terials, including crystalline and amorphous structures of metals, insulators, and semicon- 
ductors. A combination of one or more hetero junctions in a device is a heterostructure. 
As discussed in Sec. I, heterostructures having semiconductor substrates are widely used in 
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the fabrication of devices that have desired electron characteristics as functions of applied 
potentials.- 

Quantum structures may be classified as wells, wires, or dots, depending on whether the 
carriers are confined in one, two or three dimensions. Semiconductors can be used in all of 
these structures, but we consider one- dimensional potential wells and barriers, which can be 
fabricated by growing a single layer of one material (such as GaAs) between two layers of 
different materials (such as AlGaAs), as in Fig. E^iii). Other pairs of semiconductors having 
profound importance in technological applications include the III-V compounds GalnAs/InP, 
GalnAs/AlInAs, GaSb/InSb and the II- VI compounds CdZnSe/ZnSe, ZnSTeSe/ZnSSe. 

In a semiconductor the difference between the conduction and valence band energies of 
the two materials at the heterojunction with respect to the vacuum level and the Fermi level 
is responsible for the formation of quantum wells, whose heights are of the order of a few 
hundred meV. This height is to be compared with the thermal energies of carriers at room 
temperature (~ 26meV). Hence, the thermal motion of the electrons does not allow them 
to frequently cross the barriers, and quantum mechanical tunneling is the main transport 
phenomena at length scales less than the mean free path of the electrons. 

From Fig. [8p) it is evident that the conduction band potential step is given by AE C = 
Xi ~ X2- Because the Fermi level must be continuous in chemical and thermal equilibrium, 
the valence band potential step is found to be AE V = Eq2 — Eqi — AE C . 

The I-V characteristics of a one- dimensional quantum heterostructure can be related to 
the transmission probability according to the relation 2 - 

poo poo c\rp 

J= W. dk l W)-m)\T-T^, (38) 

where ki and k t are the wave vectors in the longitudinal and transverse directions respectively, 
and / is the density of states given by the Fermi-Dirac distribution. The theoretical results 
can be compared with experimental data available from the current-voltage characteristics 
of the device. The mass to be used in the Schrodinger equation is the effective mass of 
the electron m* in the longitudinal direction. As mentioned in Ref. |23| the electrons lose 



coherence after tunneling through a distance of the order of the mean free path of the 
carriers, resulting in a widening of the peaks in the I-V characteristics. Calculations of the 
bound and quasi-bound states of quantum heterostructures are available.— 

There are two major approximations in the modeling of semiconductor heterostructures 
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by rectangular potentials. The effective mass m* of the carrier changes when the electrons 
pass through a heteroj unction, but this variation is not incorporated in the analysis. Also, 
the potential profile is not steplike - the energy bands of the two materials at the het- 
erojunction change smoothly because of factors such as the inhomogeneities present at the 
interface, space-dependent composition of the compounds, and possible mechanical strain of 
the layers. The best known example is the SiC^/Si heteroj unction with a very small density 
of defects at the interface. Compound semiconductors usually have a larger concentration 
of defects and inhomogeneities, leading to a continuous barrier rather than a rectangular 
one.— >2£ The experimental data indicate a double Gaussian distribution of heights for the 
potential barriers. Another major reason for the smoothening of the heterojunction potential 
shapes is the formation of a space charge. 27 

Our results may be used to compare models with experimental data which are available 
for semiconductor heterostructures. This aspect requires more work - we need to find out the 
theoretical I-V characteristics from the transmission coefficients obtained by using Eq. (1381) . 
In principle, other functional forms may also be chosen and the I-V characteristics obtained 
and compared with experiments to determine the actual distribution of barrier heights. 
It is possible that such investigations may help design device applications where precise 
knowledge about fabricated structures is often very useful. 

VI. REMARKS 

We have carried out similar computations for the Lorentz barrier, and obtained results 
closer to the rectangular barrier than to the Gaussian barrier. Other functional forms of the 
barriers may be studied to learn their characteristic features. A further challenging problem 
is to consider multiple barriers and external, applied field effects, which are necessary to 
model realistic systems in condensed matter /semi conductor physics. 



20 



Acknowledgments 

We express our gratitude to the anonymous reviewers for their valuable suggestions in 
improving the article. 



Electronic address: quantumavik@gmail.com 



' Electronic address: sayan@iitkgp.ac.in 



1 A. Beiser, Concepts of Modern Physics (McGraw-Hill International, Singapore, 1984). See Chap. 
5 for a discussion on the tunneling effect and Chap. 12 for its application to a-decay. 

2 E. Merzbacher, Quantum Mechanics (John Wiley & Sons, New York, 1986). 

3 C. Cohen- Tannoudji, B. Liu, and F. Laloe, Quantum Mechanics, Vol. 1 (John Wiley & Sons, 
New York, 1977). 

4 V. V. Mitin, V. A. Kochelap, and M. A. Stoscio, Quantum Retero structures: Microelectronics 
and Optoelectronics (Cambridge University Press, Cambridge, UK, 1999), pp. 73-87. 

5 Z. I. Alferov, "The double heterostructure concept and its applications to physics, electronics 
and technology," Rev. Mod. Phys. 73(3), 767-782 (2001). 

6 P. Descouvemont, D. Baye, and P.-H. Heenen, "A barrier penetration model for heavy ion fusion 
valid at all energies," Z. Phys. A 306, 79-88 (1982). 

7 L. Randall and R. Sundrum, "An alternative to compactification," Phys. Rev. Lett. 83, 4690- 
4693 (1999) ibid. "A large mass hierarchy from a small extra dimension," Phys. Rev. Lett. 83, 
3370-3373 (1999). 

8 K. Rapedius and H. J. Korsch, "Barrier transmission for the one-dimensional nonlinear 
Schrodinger equation: Resonances and transmission profiles," Phys. Rev. A 77, 063610 (2008). 

9 H. A. Ishkhanyan and V. P. Krainov, "Above barrier reflection of cold atoms by resonant 
laser light within the Gross-Pitaevskii approximation," Laser Physics 19, 1729-1734 (2009); H. 
A. Ishkhanyan and V. P. Krainov, "Multiple-scale analysis for resonance reflection by a one- 
dimensional rectangular barrier in the Gross-Pitaevskii problem," Phys. Rev. A80, (045601-1-4) 
(2009); H. A. Ishkhanyan, V. P. Krainov, and A. M. Ishkhanyan, "Transmission resonances in 
above-barrier reflection of ultra-cold atoms by the Rosen-Morse potential," J. Phys. B: At. Mol. 
Opt. Phys. 43, 085306 (2010). 



21 



10 T. Paul, K. Richter, and P. Schlagheck, "Nonlinear resonant transport of Bose-Einstein con- 
densates," Phys. Rev. Lett. 94, 020404(2005). 

11 B. Mendez and F. Dominguez-Adame, "Numerical study of electron tunneling through het- 
erostructures," Am. J. Phys. 62, 143-147 (1994). 

12 A. P. Stamp and G. C. Mcintosh, "A time-dependent study of resonant tunneling through a 
double barrier," Am. J. Phys. 64, 264-276 (1996). 

13 A. U. Maheswari, P. Prema, S. Mahadevan, and C.S. Shastry, "Quasi-bound states, resonance 
tunnelling, and tunnelling times generated by twin symmetric barriers," Pramana J. Phys. 
73(6), 969-988 (2009). 

14 H. Yamamoto, "Resonant tunneling condition and transmission coefficient in a symmetrical 
one-dimensional rectangular double-barrier system," Appl. Phys. A 42, 245-248 (1987). 

15 B. Ricco, M. Ya. Azbel, "Physics of resonant tunneling: The one-dimensional double barrier 
case," Phys. Rev. B 29, 1970-1981 (1984). 

16 J. L. Powell and B. Crasemann, Quantum Mechanics (Addison- Wesley, Reading, MA, 1965). 

17 P. C. W. Davies, "Quantum tunneling time," Am. J. Phys. 73, 23-27 (2005). 

18 T. Hartman, "Tunneling of a wave packet," J. Appl. Phys. 33, 3427-3433 (1962). 

19 E. H. Hauge and J. A. Stoveng, "Tunneling times: A critical review," Rev. Mod. Phys. 61, 
917-936 (1989) 

20 V. S. Olkhovsky, E. Recami, and G. Salesi, "Superluminal tunneling through two succesive 
barriers," Europhys. Lett., 57, 879-884 (2002). 

21 V. Petrillo and V. S. Olkhovsky, "Time asymptotic expansion of the tunneled wave function for 
a double barrier potential," Europhys. Lett. 74, 327-333 (2006). 

22 H. G. Winful, "Tunneling time, the Hartman effect and superluminality: A proposed resolution 
of an old paradox," Physics Reports 436, 1-69 (2006). 

23 R. Tsu and L. Esaki, "Tunneling in a finite superlattice," Appl. Phys. Lett. 22, 562-564 (1973). 

24 E. Anemogiannis, E. N. Glytsis, amd T. K. Gaylord, "Bound and quasibound state calculations 
for biased/unbiased semiconductor quantum heterostructures," IEEE J. Quantum Electronics 
29, 2731-2740 (1993). 

25 S. Chand and J. Kumar, "Evidence for the double distribution of barrier heights in Pd,2Si/n—Si 
Schottky diodes from I-V-T measurements," Semicond. Sci. Technol. 11, 1203-1208 (1996). 

26 S. Zeyrek, M. M. Bulbul, S. Altindal, M. C. Baykul, and H. Yuzer, "The double Gaussian 

22 



distribution of inhomogeneous barrier heights in Al/GaN/p-GaAs (MIS) Schottky diodes in 
wide temperature range," Braz. J. Phys. 38(4), 591-597 (2008). 

L. V. Chebotarev, "Resonant delay in tunnelling through smooth double and triple barriers," 
Phys. Stat. Sol. (b) 208, 69-86 (1998). 



23 



